Resonances for activity waves in spherical mean field dynamos 
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We study activity waves of the kind that determine cyclic magnetic activity of various stars, 
including the Sun, as a more general physical rather than a purely astronomical problem. We try 
to identify resonances which are expected to occur when a mean-field dynamo excites waves of 
quasi-stationary magnetic field in two distinct spherical layers. We isolate some features that can be 
associated with resonances in the profiles of energy or frequency plotted versus a dynamo governing 
parameter. Rather unexpectedly however the resonances in spherical dynamos take a much less 
spectacular form than resonances in many more familiar branches of physics. In particular, we find 
that the magnitudes of resonant phenomena are much smaller than seem detectable by astronomical 
observations, and plausibly any related effects in laboratory dynamo experiments (which of course 
are not in gravitating spheres!) are also small. We discuss specific features relevant to resonant 
phenomena in spherical dynamos, and find parametric resonance to be the most pronounced type of 
resonance phenomena. Resonance conditions for these dynamo wave resonances are rather different 
from those for more conventional branches of physics. We suggest that the relative insignificance 
of the phenomenon in this case is because the phenomena of excitation and propagation of the 
activity waves are not well-separated from each other and this, together with the nonlinear nature 
of more-or-less realistic dynamos, suppress the resonances and makes them much less pronounced 
than resonant effects, for example in optics,. 



PACS numbers: 95.30.Qd, 97.10.Jb 



I. INTRODUCTION 



The well known solar activity cycle is associated with 
propagation of a wave of quasistationary magnetic field 
somewhere within the solar convective envelope. This 
wave is believed to be excited and supported by electro- 
magnetic induction effects driven by the joint action of 
differential rotation and mirror-asymmetric flows. This 
is known as dynamo action, e.g. Correspondingly, 
the wave is known as a dynamo wave. More specifically, 
there are at least two, more or less independent, dynamo 
waves, one propagating equatorwards in each solar hemi- 
sphere. 

For quite a long time, solar activity waves provided 
the unique example of dynamo wave propagation, and so 
the phenomenon was not investigated in a wider context. 
Even such a limited statement of the problem provides a 
quite rich range of instructive physical phenomena. From 
time to time the dynamo engine " misfires" , and solar 
Grand Minima, including the Maunder Minimum in the 
XVIIth - early XVIIIth centuries, occur (see for review 
e.g. [2]). At the end of the Maunder Minimum, the 
dynamo wave propagated in one solar hemisphere only, 
and a so-called "mixed-parity" magnetic configuration 
appeared [|[ . Analysis of XVIIIth century sunspot data 
[i[ gives a hint [f| that from time to time magnetic field 
becomes concentrated near to the solar equator, rather 
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than in two propagating belts in middle latitudes, and 
a configuration of quadrupole-like (instead of the usual 
dipole-like) symmetry then appears. 

Recently the variety of physical phenomenon associ- 
ated with dynamo wave propagation has enlarged sub- 
stantially. The first attempts to isolate dynamo wave 
propagation from stellar activity data were undertaken in 
[a, 0] , and generation of oscillating magnetic field, prob- 
ably indicating dynamo wave propagation, was obtained 
in the VKS dynamo experiment |8| . This opens new per- 
spectives in the topic, in particular the opportunity to 
study dynamo waves as a physical, rather than a purely 
astronomical, phenomenon. 

This change of viewpoint in the investigation of dy- 
namo waves has attracted attention to these problems, 
which had previously remained underinvestigated. In 
particular, it was stressed [9| that resonance, which is 
a prominent and sometimes spectacular phenomenon in 
the conventional theory of various oscillations and wave 
propagations, must also play an important role in dy- 
namo wave propagation. We share their view that the 
topic was not adequately addressed in previous studies, 
and only a few papers (including [l(], E|) nave addressed 
the topic. The point is that the natural statement of the 
problem relevant to astronomical contexts did not focus 
attention on this issue. 

Here we present our findings associated with the prob- 
lem of resonance phenomena in dynamo wave propaga- 
tion. Our general conclusion can be formulated as fol- 
lows. Resonances in dynamo systems do occur, and some- 
times they can be isolated in the solutions of the govern- 
ing equations (e.g. [Hi El). However in general dynamo 
wave propagation as a physical phenomenon differs quite 
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substantially from, say, the propagation of acoustic or 
electromagnetic waves. The point is that conventional 
wave propagation can be easily separated from the wave 
excitation, whereas dynamo wave propagation is usually 
intimately associated with dynamo wave excitation. As 
a result, solutions of the dynamo equations usually fail 
to separate in any explicit way resonant phenomena from 
the more general phenomena of dynamo wave excitation. 

We illustrate the problem by numerical investigation 
of a dynamo model in which dynamo waves are excited 
in two adjacent spherical shells with different distribu- 
tions of dynamo governing parameters. The magnetic 
field is continuous across the boundary between the two 
domains. A related model was discussed previously in 
fl2l [l3[ in a more specifically astronomical context. 



II. THE DYNAMO MODEL 

We consider a standard mean field dynamo equation 
(a 2 w-dynamo) 



<9B 

— = curl (VxB 
dt K 



aB — rj curl B), 



(1) 



where r\ is the turbulent diffusivity and a represents the 
usual isotropic alpha-effect (for details of the model see 
121: the codes used there are in turn based on that of 
14j). In Eq. (JlJ V is the velocity from the (differential) 
rotation and we only look for axisymmetric solutions. 
A simple algebraic a-quenching (in which a is reduced 
by a factor 1 + (B/B ) 2 , where B is a field strength at 
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which dynamo saturation is expected) is used to suppress 
dynamo action as the magnetic field grows. 

We solve the equation in spherical geometry for two 
adjacent dynamo active layers, the 'deep' layer I with 
0.7 < r < 0.85 and the 'shallow' layer II with 0.85 < r < 

I (r is the fractional radius). We consider a very simple 
rotation profile with rotation velocity linear in r, and a- 
profiles in each layer are proportional to cos 6 where 8 is 
the polar angle, so that 9 = 90° is the equator. C ! a and 
C^ 1 are dimensionless amplitudes of a in layers I and 

II correspondingly, and C<f, and C^ 1 are dimensionless 
amplitudes of the differential rotation. We take = 
C^J in the examples discussed below, a is continuous 
across the interface at r = 0.85. All solutions discussed 
have odd (" dipole-like" ) parity. 

Our investigation is focussed on the more physical as- 
pects of the problem, but it would be straightforward to 
rescale the results for particular celestial bodies, or to 
take a more realistic rotation curve. 



III. LOOKING FOR RESONANCE IN DYNAMO 
SOLUTIONS 

A. Background 

We obtained solutions for our model for various gov- 
erning parameters. We first determined the periods when 
dynamo action occurs in one layer only, i.e. either C T a or 
C^ 1 is zero. Note that the magnetic field does penetrate 
the passive layer, although its strength is substantially 
reduced there. Here we either take > so that with 
C„ = we obtain poleward migration, or C^ 1 < so 
that with C l a = there is equatorward migration. The 
variations of period with |C^' /7 | for these cases are shown 
in Fig. 1. 
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FIG. 1: Periods of oscillation in energy when dynamo action 
occurs in a single layer. Crosses: C ! a = 0, C^ 1 varies; aster- 
isks: Cl 1 = 0, Cl varies . 



FIG. 2: Time variations of magnetic energy for = 6, GjJ = 
—5. Upper panel - toroidal field energy in lower (solid curve) 
and upper (broken) layers separately, lower panel - total mag- 
netic energy. 
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More generally, the model includes two interacting dy- 
namo regions, which are physically separate but adja- 
cent, which possess specific individual eigenfrequencies, 
so quite naturally the resulting time-series for the total 
magnetic energy and magnetic energy for the upper and 
lower layer can be more complicated than just a harmonic 
oscillation; see Fig. 2 where the result for a typical solu- 
tion (Cq = 6, C„ — —5) is shown. 



B. Resonant solutions 

We then set C^ 1 — — 5 and allowed to vary, and 
plotted the period of the total energy (defined by a mean 
over 0(1000) oscillations) The variations of toroidal field 
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FIG. 3: The case C" = -5, C' u = C*" = 10 5 . Top panel: 
toroidal field energies in the lower and upper layers (asterisks 
and crosses respectively). The lower panel shows the total 
energy and the middle panel shows the variation of the pe- 
riod of oscillation of the energy. In each panel, quantities are 
plotted as a function of C7„. 



energies in the upper and lower layers, the (common) pe- 
riod and the total magnetic energy are shown in Fig. 3 
as a function of C^. We see a feature, in the form of a 
local peak of period, in the region C,£ = 5 ± 1. A related 
feature appears in the plots of all three quantities (al- 
though it is maybe less pronounced in the lower panel). 
Examination of Fig. 2 shows that this feature does not 
correspond to values for which the periods of the lay- 
ers, taken individually, are equal. For C^ 1 — —5, periods 
are approximately equal when « 7. 

We show in Fig. 4 short samples of the time series 
of the total energy for solutions outside of the resonant 
domain (C^ = 6.5, upper panel), and inside the domaon 
(C„ = 4.5, lower panel). 



C. Resonance with a subcritical dynamo layer 

We recognize that the resonant behaviour was not very 
well pronounced in the previous results and thus try to 
explore and explain this situation. We do not doubt 
that resonance in its purely mathematical sense occurs 
in dynamo solutions. More specifically, if the govern- 
ing equations for a kinematic dynamo problem have a 
degenerate eigenvalue r n = 7„ + ioJ n (7n > 0), then 
it is more than plausible that the corresponding set of 
eigenvectors becomes insufficient to provide a basis, and 
the desired solution, apart from terms proportional to 
exp(7„ + iu> n )t, may contain resonant terms proportional 
to iexp(7„ + iu) n )t. The point however is that the ex- 
ponential growth of the solution as exp7„i or texpj n t 
occurs just at the initial stage of the dynamo wave ex- 




FIG. 4: Extracts from the time series of total energy for solu- 
tions outside of the resonant domain (upper panel, C ! a = 6.5), 
and inside this domain (lower panel, Ci = 4.5). 
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citation and is a matter of more or less purely theoret- 
ical interest. A more practical issue is the amplitude 
of the steady-state oscillation obtained after the expo- 
nential growth is somehow saturated - here by the term 
1 + (B/Bq) 2 . The saturation reduces the growth rate 
7„ to zero, and this is a much stronger effect than the 
additional power-low growth due to the resonance. The 
exceptional situation discussed by, say, |10Ulll |) addresses 
the cases where excitation of a particular solution with- 
out a resonance does not occur, and it is the existence 
of a resonance that allows this solution to develop. In 
contrast, resonances in acoustic or electromagnetic wave 
propagation occur when a wave can marginally propagate 
and does not grow (7 = 0), and the resonance is clearly 
recognisable as a growing (and then saturated) solution. 

In another investigation, the lower layer was made 
much thinner (0.7 < r < 0.775), so a dynamo in this 
layer alone is subcritical (i.e. with C ! a = < 10, C^ 1 = 0, 



dynamo action does not occur) and the upper layer corre- 



spondingly thicker, 0.775 < r < 1.0. We keep C a 
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FIG. 5: The case C' a = -5, d = C" = 10 5 , with thin lower 
layer, 0.7 < r < 0.775. Top panel: toroidal field energies in 
the lower and upper layers (asterisks and crosses respectively). 
The lower panel shows the total energy and the middle panel 
shows the variation of period. 



C^, = C*J = 10 5 . There is again a (smaller) feature 
attributable to a resonance near C J a — 4.25. See Fig. 5. 

We conclude from Fig. 5 that the resonance here is of 
a similar nature to that in the previous case. 



D. Parametric resonance 

Now we seek the other kind of resonance, i.e. a para- 
metric resonance in the dynamo wave driven by harmonic 
variations of the dynamo governing parameters. Our pre- 
liminary expectations here are twofold. From one hand, 
the resonances recognized in [l(| [DJ for galactic dynamos 
are parametric. On the other hand, the parametric res- 
onance suggested in [l5[ to explain some details in the 
magnetic activity of stellar binary systems does not seem 
very spectacular. 

We imposed a modulation a = a (l + fcosujpt) in 
the region I, with 500 < lo p < 1500, / = 0.2, keep- 
ing C^ 1 — — 5. We started from the unperturbed model 
which has C J a — 7.5 (i.e. outside of the resonant re- 
gion discussed in Sect. IIII B|) . The corresponding period 
is P{E) = 0.00539, so w{E) = 1166. Results are shown 
in Fig. 6, where resonant peaks are clearly recognizable. 

The feature at ujp « 1150 is quite expected, as it sat- 
isfies the standard resonance condition ujp k, oj{E) — 
2lj(B) (where oj(B) is the frequency of the magnetic field 
oscillations). The features near lu p 700 and oj p w 900 
have to be associated with rather unusual excitation con- 
ditions. Note that these resonant solutions do not settle 
to a completely steady (maybe modulated) oscillation, 
even when run for about 2000 oscillations. In the two 
anomalous cases the period is not quite stable, but the av- 
eraging process seems robust, in that averages over differ- 
ent temporal sub-ranges give the same values. Note that 
the energy in the upper and lower regions separately (top 
panel in Fig. 6) is the toroidal field energy only, whereas 
the "total energy" includes also the poloidal. Over long 
enough time intervals the average period in region II is 
slightly longer than that of the total energy (see middle 
panel of Fig. 6). 

In general, parametric resonance effects seem to be 
much clearer and more pronounced than the resonances 
analysed in previous models fSects. lTlI Bl and llll Cp . How- 
ever the peaks in energy profiles for parametric reso- 
nances are quite broad. 



IV. DISCUSSION 

The concept of resonance belongs to the basic concepts 
of contemporary physics, and it looks a priori implau- 
sible that it should not be applicable to dynamo wave 
propagation. We have verified this natural expectation 
in the framework of activity waves generated by dynamo 
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waves, in a quite traditional model of a two-layer mean- 
field spherical mean-field dynamo of the type convention- 
ally considered as a model for periodic stellar magnetic 
activity. 

Indeed, we found features in the profiles for energies 
or periods (frequencies) plotted versus a dynamo govern- 
ing parameter which can be identified with a resonance. 
However, it looks noteworthy that the details resemble 
more discontinuities in the profiles of the behaviour of 
period and energy as varies, rather than the con- 
ventional resonant peaks shown in physics textbooks. It 
is even difficult to claim that the resonance makes dy- 
namo excitation more efficient and enhances, say, the 
total magnetic energy. The lower panel in Fig. 5 is in- 
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FIG. 6: Parametric resonance: notation as in Fig. 3. The 
middle panel now includes the periods of oscillations in the 
upper layer II (shown by small solid circles). The dashed 
horizontal lines are values for the unperturbed model. The 
period of the total energy is plotted by solid dots in the middle 
panel. In each panel, quantities are plotted a functions of 
W (P). 



structive in this respect: the interaction here appears as 
a fall in total energy, rather than an increase. Of course, 
resonances are far from being the only effects that de- 
termine the relevant behaviour in, say, optics, and there 
are many cases where other effects make recognition of 
resonance peaks problematic. The point however is that 
it is not difficult to provide examples of pronounced res- 
onant peaks in optics. When presenting our results, we 
have of course chosen the more pronounced examples of 
resonances for dynamo waves, admittedly from a rather 
small sample. 

Correspondingly, we cannot suggest that the sort of 
resonances discussed here will leave observable signa- 
tures. Our intention is to draw attention to novel be- 
haviour. Our results are not so marked as those in 
arguably because a linear problem was considered there. 
The inclusion of a nonlinear saturation mechanism re- 
duces these effects. 

The parametric resonances present in the spherical dy- 
namo models investigated here (Sect. IIII D[) are much 
clearer than the other resonance features (Sects. IIII Bl 
and IIII Cj) . Perhaps, this is a consequence of controlling 
by varying the frequency u)p, which governs the para- 
metric excitation but, leaves unchanged the "nominal" 
parameters of the dynamo model for the vanishing mod- 
ulation amplitude / = 0. 

Our numerical experiment illuminate the unusual ex- 
citation conditions for parametric resonances. Perhaps, 
it can be said that the Mathieu equation is not a fully 
adequate model for parametric resonance in spherical dy- 
namos. Isolation of the resonance condition here appears 
an attractive topic for dynamo theory, but it is however 
obviously beyond the scope of this paper. The transi- 
tion time required for the dynamo system to achieve a 
resonant type solution can be long and a qualitative ex- 
planation in terms of an appropriate model equation also 
looks desirable. 

We note also the existence of a quasi-resonant phe- 
nomenon in a case where one of the dynamo layers is 
passive, i.e. with the chosen parameters the lower layer 
taken in solution does not support dynamo action. 

Perhaps the basic differences between resonant phe- 
nomena for dynamo waves and those in more conven- 
tional branches of physics can be summarized as fol- 
lows. It is difficult to separate dynamo wave propagation 
from dynamo wave excitation, while such processes are 
well separated in many other branches of physics. Cor- 
respondingly, eigenvalues of the linear (kinematic) dy- 
namo problems take complex values. Coincidence of two 
complex quantities is a more severe requirement rather 
than equality of two real valued quantities. Dynamo self- 
excitation is in practice saturated by some back-reaction 
of the dynamo generated magnetic field on the hydrody- 
namics. This saturation determines the resulting mag- 
netic energy and affects the periods of dynamo waves. 
The corresponding effects specific to dynamo action are 
more important than any potential resonant behaviour. 
In additional, the two layers of dynamo action are sep- 
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arated hydrodynamically but magnetic field penetrates 
from one activity layer to the other, even if the dynamo 
action in one layer remains subcritical. Thus the very 
concept of the frequency of dynamo waves in one activ- 
ity layer taken alone is not fully applicable. 

In general, it seems that the situation under discus- 
sion reflects the every-day wisdom that small problems 
(here resonances) can only spoil your life in the absence 
of larger ones (in this case other uncertainties associated 



with mean field dynamos are much larger). 
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